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Abstract 

In this note I introduce a mysterious approximation called the rotating wave approx- 
imation (RWA) to undergraduates or non-experts who are interested in both Mathe- 
matics and Quantum Optics. 

In Quantum Optics it plays a very important role to obtain an analytic approximate 
solution of some Schrodinger equation, while it is curious from the mathematical point 
of view. 

I explain it carefully with two coherent oscillations for them and expect that they 
will overcome the problem in the near future. 
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1 Introduction 



When undergraduates study Quantum Mechanics they encounter several approximation 
methods hke the WKB, the Born-Oppenheimer, the Hartree-Fock, etc. In fact, exactly 
solvable models are very few in Quantum Mechanics, so (many) approximation methods 
play an important role. As a text book of QM I recommend [1] although it is not necessarily 
standard. 

When we study Quantum Optics we again encounter the same situation. We often 
encounter a method called the rotating wave approximation (RWA), which means fast oscil- 
lating terms (in effective Hamiltonians) removed. Because 



J ±in 

if n is large enough. We believe that there is no problem. 

However, in some models slow oscillating terms are removed. Let us show an example. 
The Euler formula gives 



because e~*^ goes away from e*^ by two times speed, so we neglect this term. In our case n 
is 2 ! Read the text for more details. 

Why is such a "rude" method used ? The main reason is to obtain analytic approximate 
solutions for some important models in Quantum Optics. To the best of our knowledge we 
cannot obtain such solutions without RWA. 

In this review note I introduce the rotating wave approximation in details with two 
models for undergraduates or non-experts. I expect that they will overcome this "high wall" 
in the near future. 





From this we approximate 2 cos 6 to be 



2 cos ^ = e*^ + e 
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2 Mathematical Preliminaries 



We make a short review of the two-dimensional complex vector space and complex matrix 
space M(2; C) within our necessity in order to treat the two-level system of an atom. See 
for example [2]. 

First we introduce the (famous) Pauli matrices {cri, (72, cts} defined by 



1 

1 

and set the unit matrx I2 by 



0-2 



-i 

1 

1 
1 



era 



1 
-1 



Moreover, we set 



1 




0"_ 





1 



Note that (Xi = (T_|_ + (T_ . Then it is easy to see 



[^0-3,0-+] = 0-+, [io-3,cr_] = -cr_, [o-+,0-_] = 2 X i(T3. 



(2) 



Comment The Pauli matrices {cti, (T2, (J3} are generators of the Lie algebra su{2) of the 
special unitary group SU{2) and {(j+,cr_, ^cr^} are generators of the Lie algebra sZ(2; C) of 
the special linear group SL{2] C). For the sake of readers we write a Lie-diagram of these 
algebras and groups. 



s/(2;C) 



su{2) 



SL{2;C] 



SU{2) 



Next, we define {|0), |1)} a basis of by use of the Dirac's notation 



|0) 



, |1) 





1 



(3) 



Then, since di satisfies the relation 

ai|0) = |l), ai|l) = |0) 

it is called the flip operation. 

Note If we deflne {(J+, cr^, ^0-3} as above then {|0), |1)} should be chosen as 
instead of ([3j). Because, 

a_|0) = 0, fT_|l) = |0), a+|0) = |l). 
However, I use the conventional notations in this note. 

For the later convenience we calculate the exponential map. For a square matrix A the 
exponential map is deflned by 



00 



n\ ^ — ' n! 

ra=0 n=0 



where E is the unit matrix and A is a constant. 

Here, let us calculate e*'^'^^ as an example. Noting 



1 \ 1 

10/ \ 1 



we obtain 



n=0 



(2n)\ ^(2n + l,. 

n=0 ^ ' n=0 ^ ' 

~ (_1)"A2" (_i)-A2"+i 

n=0 ^ ^ n=0 ^ ^ 

COS A I2 + z sin A ai 
cos A i sin A 
i sin A cos A 



(4) 



Exercise Calculate 



e'^"^ and e'^''^ 



3 Two— Level System of an Atom 



We treat an atom trapped in a cavity and consider only two energy states, namely (in 
our case) the ground state and first excited state. That is, all the remaining states are 
neglected. This is usually called the two-level approximation. See for example [3] as a 
general introduction. 

We set that energies of the ground state |0) and first excited state |1) are Eq and EiiEq < 
E^ijrespectively. Under this approximation the space of all states is two-dimensional, so there 
is no problem to identify {|0), |1)} with ([3]). 

Then we can write the Hamiltonian in a diagonal form like 

Eq 

El 

For the later convenience let us transform it. For A = Ei — Eq the energy difference we have 



Ho 



Eq+Ei _ E^-Eq 
2 2 



2 2 



Eq + Ei ^ _ El — Eq ^ 
Eo + Ei A 



(5) 



To this atom we subject LASER (Light Amplification by Stimulated Emission of Radiation) 
in order to controU it. As an image see the following figure. 

El |1) ^ ^ 

|0) J J 



In this note we treat L cLSGr clS db classical wave for simplicity, which is not so bad as shown 
in the following. That is, we may set the laser field as 

A cos (cut + (p). 

By the way, from several experiments we know that an atom subjected by Laser raises 
an energy level and vice versa. This is expressed by the property of the Pauli matrix cti 

ai|0) = |l), ai|l) = |0), 

so we can use cxi as the interaction term of the Hamiltonian. 

As a result our Hamiltonian (effective Hamiltonian) can be written as 

H = ——CTs + 2g cos{ujt + (f))(Ji 

— Y 2(7 cos (cut + 0) 

2g cos{ujt + 0) ^ 



(6) 



where 5^ is a coupling constant regarding an interaction of between an atom and laser, and 
A is absorbed in g {gA — > g). We ignore the scalar term ^o+-^i for simplicity. Note that 
dS]) is semi-classical and time-dependent. 

Therefore, our task is to solve the Schrodinger equation 



d 

ih—m = (7) 

dt ^ ' 



exactly (if possible) 



4 Rotating Wave Approximation 

Unfortunately we cannot solve ([7]) exactly at the present time. It must be non-integrable 
although we don't know the proof (see the appendix). Therefore we must apply some ap- 
proximate method in order to obtain an analytic approximate solution. Now we explain a 
method called the Rotating Wave Approximation (RWA). Let us recall the Euler formula 

e'^ = cos^ + isin^ =^ 2 cos^ = e^^ + e'^^ 
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From this we approximate 2 cos 9 to be 



-2ie\ ^ je 



because e *^ goes away from e*^ by two times speed, so we neglect this term ! We call this 
the rotating wave approximation. 




Problem In general, fast oscillating terms may be neglected because 







if n is large. Our question is : Is n = 2 large enough ? 



By noting that the Hamiltonian is hermitian, we approximate 



2 COs{0Jt + 0)(Ti 



2 cos(a;t + 0) 

2 cos(a;t + 0) 

by use of ([8]), so ([6]) becomes 



H 



2 



As a result our modified task is to solve the Schrodinger equation 

d 

ih—^ = Hm 
dt 

exactly. Mysteriously enough, this equation can be solved easily. 



-i(Lut+<j>) 



(9) 



(10) 



Note For the latter convenience let us rewrite the method with formal notations : 



2 cos 61 (Ti = cTi (g) 2cos6' = ((7+ + cr_) (g) (e*^ + e 



0-+ (g) e^^ + a+ (g) e + a_ (g e*^ + cr_ (g e"*^ — ^ ct+ (g e*^ + o-_ (g e~^^ 



In order to solve (ITU]) we set ^ = 1 for simplicity. From ([H]) it is easy to see 
-f ^e^(-+^) \ _ / W -f ^ W 

so we transforn the wave function in ffTOl) into 



$=|^ ' |^^^=(^' 1$. (11) 



Then the Schrodinger equation f[TU|) becomes 



^^$=f ''''' ^ 1$ (12) 



2 



by a straightforward calculation. 

Here we set the resonance condition 

A = oj (<^=^ El — Eq = hu precisely). (13) 

Namely, we subject the laser field with u equal to the energy difference A. See the following 
figure. 



U! El 



El — Eq = hoj 

|0) 



Then (IT^ becomes 



dt 



d g 

i— $ =1 I $ = gai^ 

9 
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and we have only to solve the equation 

d 



= -tg(Xi<l>. 



By dll) (A = —gt) the solution is 



cos(qt) —i sm(qt) , 
= e-*f*"i<l>(0) = I I $(0), 

—ism{gt) cos{gt) 

and coming back to (from $) we obtain 

e*^^^ \ / cos(qt) —ism(qt) 



e * 2 y y — ■isin(5ft) cos{gt) 
g-»M+<^) / \ _i;sin(5(t) cos{gt) 



cos(5ft) — 2sin((yft) 

sm{gt) e-*('^*+'^) cos(^t) 



M/(0) (14) 



by f lTT|) (^l'(O) = $(0)) because the total phase e** 2'*' can be neglected in Quantum Me- 
chanics. 

As an initial condition, if we choose 



^(0) = |0) 







we have 



cos{gt) \ _ I cos{gt) 
_^^-i{ut+<t>) sm{gt) I \ e-*M+'^+^/2) sin(^t) 



cos{gt)\0) + e-*(-*+^+-/2) sm{gt) \ 1) (15) 



by (jSj). That is, \E'(t) oscillates between the two states |0) and |1). This is called the coherent 
oscillation or the Rabi oscillation, which plays an essential role in Quantum Optics. 

Concerning an application of this oscillation to Quantum Computation see for example 

El- 
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Problem Our real target is to solve the Schrodinger equation 

d 
dt 

with 

(— 2qcos(ut + 
2 
2g cos{ut + (p) f 

Present a new idea and solve the equation. 



5 Quantum Rabi Model 

In this section we treat the quantum Rabi model whose Hamiltonian is given by 

H = — (T3 ® 1 + UI2 ® a^a + gai (a + a') 

= ^a3®l+ul20N + g{a+ + a_)®{a + a^) (16) 

where 1 is the identity operator on the Fock space J-" generated by the Heisenberg algebra 
{a, a"'',A^ = a^a}, and Q and u are constant, and g is a. coupling constant. As a general 
introduction to this model see for example see |3|. 

Let us recall the fundamental relations of the Heisenberg algebra 

[N, a^] = a\ [N, a] = -a, [a, a^] = 1. (17) 

Here, the Fock space J-" is a Hilbert space over C given by 

J- = Vectc{|0),|l),--- ,|n),---} 

where |0) is the vacuum (a|0) = 0) and 

fat)" 

\n) = ^j^\0) for n>0. 
vn! 
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On this space the operators (=infinite dimensional matrices) a) , a and are represented as 



1 

y/2 

73 


/n 



N = ah 



■ J 



by use of f|T7|) . 



Note We can add a phase to {a, a"!"} hke 





1 
V2 

Vs 



b = e^^a, = e-'^a\ N = b^b = ah 



\ 



■ J 



:i8) 



where 6 is constant. Then we have another Heisenberg algebra 

[N,b^] = b\ [N,b] = -b, M = l. 
Again, we would like to solve Schrodinger equation = 1 for simplicity) 

i^l^) = i/|^) = l^cTs^l+uh^N + g{a+ + a_)®{a + a^)\\^) 



(19) 



exactly. To the best of our knowledge the exact solution has not been known, so we must 
make some approximation in order to obatin an analytic solution. 
Since 

(cr+ + a_) (g) (a + a"!") = ® a + ® a^ + a- ® a + a- ^ a\ 
we neglect the middle terms cr+ (g) a"!" + cr_ ® a and set 



H = ^0-3 (g) 1 + uh (S) N + g{a+ (g) a + (g) a^). 



(20) 
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^In [3] it is called the Jaynes-Cummings-Paul one 
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(21) 



This is called the rotating wave approximation and the resultant Hamiltonian is called the 
Jaynes-Cummings oneQ, [5]. 

Therefore, our modified task is to solve the Schrodinger equation 

exactly. Mysteriously enough, to solve the equation is very easy. 
For a unitary operator U = U{t) we set 

1$) = 

Then it is easy to see 
by (121]) . If we choose U as 



(we use Y ill place of ^1 for simplicity), a straightforward calculation gives 

UHU-' + ^^U-'=( I (22) 

and we have a simple equation 

Note that in the process of calculation we have used the relations 

The proof is easy by use of the formula 

e^Ae-^ = A+[X,A] + [X, A]] + [X, [X, A]]] + ■ • • (24) 



for square matrices X, A and ( IT7I) . 
Here we set the resonance condition 



Q = u, (25) 

then fl23l) becomes 

Let us solve this equation. By setting 



A 



a 



a 

t 



we calculate the term e Noting 





H 







we have 



E 



+ 1 

A'=\ 1 = 1 I (^[a,at] = l) 



-igtYA"" 



n=0 



~ / (iV + 1)" \ - (_l)n((^t)2"+l / (iV+l)"a 



n=0 



(2^)' I t^o (2n + l)! \ ^.^t 



cos(v/iVn(7t) \ _ / ^ sin(v/iVn(7t)i 

cos(ViV^t) y y ^sin(yiV^t)at 

cos(ViV+ l^ft) -i-^^ sin( ViV + lgt)a 
-i-^?,m{^/Ngt)a) cos{^/N gt) 

Therefore, the solution is given by 

, cos{^/N +lgt) -i^^sm{^/N +lgt)a . 

Mi)) =1 I l-^W) 

-2^sin(A/iV5(t)at cos{y/N gt) 

13 



(26) 



and coming back to |\E') (from |$)) we finally obtain 

. e**(^^+t) \ / cos(^/P^TTgt) -i^^sm(^/WTlgt)a , 

\^{t)) =1 . . J ' vWTT ^ y ^ I 1^(0)^ 

gif(^7v__) J y _^_^sin(yiV^t)at cos{VNgt) 



where |^(0)) = |$(0)). 

As an initial condition, if we choose 



we have 



— ze**2 sin(5fi)|l) 



because a|0) = and a^O) = |1), or 



|^(t)) = cos(^t) I M ® |0) + e-'^ sm{gt) P | ® |1) 



where the total phase e**" has been removed. 



Problem Our real target is to solve the Schrodinger equation 



d 



with 



H = Hit) = ^(73 (S}l+Lul2^ N + g{a+ + (7_) (g) {a + 



(27) 



l*(0)) = I J ® = ( 



. e'*^ cos(gt)\Q) . 
\^{t))= ' " ' (2^ 



Present a new idea and solve the equation. 

As a developed version of the Jaynes-Cummings model see for example |B] and [7]. 

6 Concluding Remarks 

In this note I introduced the rotating wave approximation which plays an important role 
in Quantum Optics with two examples. The problem is that the method is used even in 
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a subtle case. As far as I know it is very hard to obtain an analytic approximate solution 
without RWA. 

I don't know the reason why it is so. However, such a "temporary" method must be 
overcome in the near future. I expect that young researchers will attack and overcome this 
problem. 



Concerning a recent criticism to RWA see [8] and its references, and concerning recent 
apphcations to the dynamical Casimir effect see [9], [10] and [H], [12] ([TT] and [I2] are 
highly recommended). 



[A] Another Approach 

Let us give another approach to the derivation which may be smart enough. It is easy 
to see the diagonal form 



Appendix 



cTi = WasW 



-1 



where W is the Walsh-Hadamard matrix (operation) given by 




G 0(2). 



Note that 



= I2 



w 



w 



-1 
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Then we obtain 




Readers should remark that the Walsh-Hadamard matrix W plays an essential role in 
Quantum Computation. See for example [13] (note : — )■ [/a in this paper). 



[B] Tensor Product 

Let us give a brief introduction to the tensor product of matrices. For A = (aij) G M{m; C) 
and B = (bij) G M{n] C) the tensor product is defined by 

A(S) B = (aij) ® B = (aijB) G M{mn; C). 

Precisely, in case of m = 2 and n = 3 



A^B = 



/ o-ii ai2 \ 




Q-\\B 


ayiB \ 








\®B=\ 










\ 0-21 O22 j 




O21-B 


a-22B J 






^ aiifoii 


011612 


011613 


012611 




012613 


\ 


011621 


6^11622 


011623 


O12621 


O12622 


012623 




011631 


011632 


011633 


012631 


012632 


012633 




021611 


021612 


021613 


022611 


022612 


022613 




'^21621 


O21622 


021623 


O22621 


O22622 


O22623 




\ ^21631 


O21632 


O21633 


O22631 


O22632 


022633 


/ 



When I was a young student in Japan this product was called the Kronecker one. Nowa- 
days, it is called the tensor product in a unified manner, which may be better. 
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Note that 



while 



( 



\ 



bn 


bi2 


bi3 








&21 


&22 


&23 








&31 


^32 


&33 














7 

&11 


7 

O12 


7 

Ol3 








7 

O21 


7 

O22 


7 

O23 








7 

O31 


7 

O32 


7 

O33 


bu 




bl2 




&13 






Oil 




h 




013 


&21 




&22 




&23 










b22 




^23 


631 




&32 




&33 






bsi 




bz2 




^33 



'33 / 



The blanks in the matrices above are of course zero. 

Readers should recognize the difference. See for example [2] for more details. 



[C] Beyond the RWA 

Let us try to solve the equation ([7]). For the purpose it is convenient to assume a form for 
some solution 

where we set = (l/2)cr3 for simplicity. Note that this form called the disentangling form 
(a kind of Gauss decomposition of some matrices) is very popular in Quantum Physics. 
By setting ^ = 1 for simplicity in ([7]) we must calculate 

d 

i-^^ = {-Ar3 + 2^cos(wt + 0)(cr+ + (T_)}^ 

= {2g cos{ut + 0)(T_|_ — At3 + 2g cos{ut + 0)(T_ } \& 

where T3 = (1/2)(T3. 
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Then we have 



= ^^|g-iF(tWg-iG(t)r3g-iH(tV_|^(Q) 



dt 



From (ED 



[rg, cr+] = (7+, [ra, (7_] = -a_, [0-+, (t_] = 2r3 



and the formula fl24j) it is easy to see 



;i + zG(t))a_, e-*^w-+^_e 



Therefore 



T 



F{t)a+ + G(t)(r3 + iF{t)a+) + + iG{t)){a^ - 2iF{t)n + F(t)V+)} ^ 

+ ^F(t)G'(t) + (1 + iG{t))F{tfH{t)^ a++ 
G(t) - 2i(l + iG{t))F{t)H{t)^ ^3 + (1 + iG{t))H{t)a_ ^. 

By comparing two equations above we obtain a system of differential equations 

F{t) + iF{t)G{t) + (1 + iG{t))F{tfH{t) = 2g cos{ujt + 0), 
G{t) - 2z(l + iG it)) Fit) Hit) = -A, 
{l + iG{t))Hit) = 2gcos{ujt + <j)). 

By deforming them we have 

F{t) - iAF{t) - 2g cos(wt + = 2g cos{ut + 0), 

G{t) - Aig cos(wt + = - A, 

(1 + iG(t))if(t) = 2c/cos(wt + 0). 

This is a simple exercise for young students. 

If we can solve the first equation then we obtain solutions like 

F{t) =^ G{t) =^ Hit). 
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The first equation 

F{t) - 2g cos{ut + 0) - iAF{t) - 2g cos{ut + (j))F{tf = 

is a (famous) Riccati equation of general type. Unfortunately, we don't know how to solve 
it explicitly at the present time. 

[D] Full Calculation 

Let us give the full calculation to the equation fl23|) . We set 



B 



ga 

ga^ 



and calculate e "^^^ without assuming Q = u in (1251) . Again, noting 



{^Y + g 



{^Y + g'N + g^ 



{^Y + g^n^ 



a'a 



i^Y + g'N 



aa) = a'a + 1 = + 1) we obtain 



-UB 



exp —it 
( 



ga 



ga^ 



cos 



\ 



2 _ j^sin£\/^+g2 
2 

• sin t.nB f 



-ig- 



■a' 



• sint^^ \ 

-%g^—a ^ 
costy^+f^ j 



where we have set 



6 = Q-uj, Lp="—+g'^N 



for simplicity. See (12^ . This is a good exercise for young students. 
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